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ABSTRACT 

The loss of angular momentum due to unstable r-modes in hot young neutron stars 
has been proposed as a mechanism for achieving the spin rates inferred for young 
pulsars. One factor that could have a significant effect on the action of the r-modc 
instability is fallback of supernova remnant material. The associated accretion torque 
could potentially counteract any gravitational-wave induced spin-down, and accretion 
heating could affect the viscous damping rates and hence the instability. We discuss 
the effects of various external agents on the r-mode instability scenario within a simple 
model of supernova fallback on to a hot young magnetized neutron star. We find that 
the outcome depends strongly on the strength of the star's magnetic field. Our model 
is capable of generating spin rates for young neutron stars that accord well with initial 
spin rates inferred from pulsar observations. The combined action of r-mode instability 
and fallback appears to cause the spin rates of neutron stars born with very different 
spin rates to converge, on a timescale of about a year. The results suggest that stars 
with magnetic fields < 10 13 G could emit a detectable gravitational wave signal for 
perhaps several years after the supernova event. Stars with higher fields (magnetars) 
are unlikely to emit a detectable gravitational wave signal via the r-mode instability. 
The model also suggests that the r-mode instability could be extremely effective in 
preventing young neutron stars from going dynamically unstable to the bar-mode. 

Key words: Accretion, accretion discs - Gravitational waves - Magnetic fields - Stars: 
neutron 



MOTIVATION 



Neutron stars may, by virtue of their oscillations, act as 
sources of detectable gravitational radiation. Of particular 
interest are the r-mode oscillations, quasi-toroidal currents 
that arise in rotating stars due to the action of the Corio- 
lis force. Andersson (1998) and Friedman & Morsink (1998) 
first showed that the r-modes of a perfect fluid star suc- 
cumb to the gravitational radiation driven Chandrasekhar- 
Friedman-Schutz instability for all rates of stellar rotation. 
Subsequent results suggest that the instability might grow 
on short timescales, with positive feedback increasing the 
gravitational wave signal to potentially detectable levels 
(Owen et al 1998). The r-mode instability is likely to be 
most important in hot young rapidly rotating neutron stars. 
For such stars there appears to be a window of tempera- 
tures and spin rates in which the r-modes can grow faster 
due to the emission of gravitational waves than they are 
damped by viscosity. The consequent loss of angular mo- 
mentum provides a plausible mechanism for achieving the 
spin rates inferred for young pulsars (Andersson, Kokkotas 
& Schutz 1999). In fact, the excitement regarding the r-mode 
instability stems to a large extent from the fact that the first 
studies suggested that it would spin a newly born neutron 



star to roughly 20 ms, in good agreement with the inferred 
initial spin period of the Crab Pulsar. 

Recent efforts to improve our understanding of the un- 
stable r-modes have focused on the interior physics, aiming 
to shed light on issues such as the role of general relativ- 
ity (Lockitch, Andersson & Friedman 2001), superfluidity 
(Lindblom & Mendell 2000), the crust/core interface (Bild- 
sten & Ushomirsky 2000; Lindblom, Owen & Ushomirsky 
2000), and possible exotic states of matter like hyperons 
(Jones 2001; Lindblom & Owen 2002). As we will show in 
this paper, another factor that may significantly affect the 
r-mode instability is accretion. Neutron stars are born in a 
messy environment, surrounded by debris from the super- 
nova progenitor. In the first few months some of this mate- 
rial will fall back on to the young star. Intuitively, one would 
expect this to affect the action of the r-mode in two ways. 
Firstly, accretion torques could counteract the gravitational 
wave induced spin-down, and perhaps cause the star to spin 
up despite the presence of a large amplitude r-mode. One 
can readily estimate that accretion of a total mass AM can 
spin a neutron star up to a rotation period 



AM 
0.1 M P 
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This estimate indicates that rapid supernova fallback may 
play a key role in determining the spin evolution of a nascent 
neutron star. Secondly, rapid accretion would be associated 
with strong heating, which would affect the strongly tem- 
perature dependent viscosities, and hence the longevity and 
strength of the r-mode instability and any associated grav- 
itational wave signal. Accretion is also expected to delay 
the formation of a crust until the star has cooled to a bulk 
temperature ~ 10 9 K (Haensel 1997). This should be com- 
pared to the case of isolated neutron stars in which the crust 
is expected to form at ~ 10 10 K. This delayed crust for- 
mation could be significant since the presence of a crust 
would suppress considerably any unstable r-mode (Bildsten 
& Ushomirsky 2000; Andersson & Kokkotas 2001). A delay 
in crust formation could therefore increase the longevity of 
any gravitational wave signal. 

Our aim with this paper is to address two key questions. 
Firstly, can the combined action of supernova fallback and 
r-mode instability give rise to spin rates that match those 
inferred for young pulsars? Secondly, what are the impli- 
cations for the gravitational wave signal from the r-mode 
instability? In order to address these questions we develop 
a simple model of supernova fallback on to a hot young 
magnetized neutron star, modelling both spin evolution and 
gravitational wave emission. 



2 THE MODEL 
2.1 Spin evolution 

In order to be able to compare our results to previous work 
and assess directly the role of the induced accretion torque, 
we employ a simple phenomenological spin evolution model 
analogous to that devised by Owen et al (1998) (see also Ho 
& Lai 2000; Andersson & Kokkotas 2001). 

We assume that we can model the star's angular mo- 
mentum as the sum of bulk angular momentum J and the 
canonical angular momentum of the r-mode J c . The total 
torque on the star J can therefore be written as 



j = m + m + j c 



(i) 



where the dots indicate time-derivatives. Q is the star's an- 
gular velocity. For simplicity we will assume that the star 
does not rotate differentially despite the presence of a large 
amplitude r-mode and a strong accretion torque. I = 1MB? 
represents the moment of inertia of the star. We will model 
the star as an n — 1 polytrope, which means that I — 0.261. 
For this particular stellar model, accretion will affect the 
mass M of the star but not the radius R. 

In our analysis we will only account for radiation from 
the (dominant) I = m = 2 current multipole of the r-mode. 
The canonical angular momentum of the mode is 



Jc 



3Sla 2 JMR 2 



(2) 



where a is the mode amplitude (as defined by Owen et al 
1998) and J = 1.635 x 10~ 2 for an n = 1 polytrope. 

The r-mode is driven by gravitational radiation and 
damped by viscosity. We thus assume that the canonical 
angular momentum, which is proportional to the square of 
the perturbation, evolves according to 



Jc 1 

2J C ~ ~ V 
where 

T _1 = — I To 



-1 i — 1 i -1 



(3) 



(4) 



and Tg,b,s are the timescales associated with growth of the 
linear perturbation due to gravitational radiation back re- 
action and dissipation of energy due to bulk and shear vis- 
cosity, respectively. From equation (^) it can be seen that 
the mode is unstable (in an isolated star) if r < 0. In this 
case a small perturbation will lead to an unbounded growth 
of the mode. 

Equations ffl), (^) and ^ can be combined to give equa- 
tions for the evolution of a and fi: 



2ft + 2M 



and 

ft = 



Mft 
~M~ 



3Ja 2 ft 



It 



(5) 



(6) 



We initiate our evolutions by assuming that a small ampli- 
tude r-mode (with a = 10 -6 ) is present. Variation of the 
initial amplitude has no significant effect on the subsequent 
evolution. 

Equations (JE]) and (^|) only apply when a is small. 
To model the behaviour once the r-mode reaches the non- 
linear regime we assume that it saturates at some value a s . 
The larger the value of a 3 , the greater the fraction of the 
star's angular momentum the r-mode can carry (we consider 
0.01 < a s < 1 in this paper). When the mode is saturated 
the spin evolution is described by 



ft 



1 



3Ja„ 



21 



Mft 



(7) 



As discussed by Owen et al (1998), the bulk of gravitational 
wave emission and associated spin-down will occur during 
the saturated phase. 

Eventually the mode will become stable again, de- 
saturate and die away; when this happens we switch back 
to using equations (|E]) and ([]). Note that equation (^) indi- 
cates that the r-mode is unstable, in the sense that a > 0, 
whenever 



M 



~ + ^ + < 0. 



(8) 



i ft 

r 2ft ' 2M 

This is the appropriate criterion for the onset of instability 
in a star that is spun up by accretion (Ho & Lai 2000). 
Compare this with the instability criterion for an isolated 
star, r < 0. 

We use the estimates given by Andersson & Kokkotas 
(2001) for the gravitational radiation reaction timescale due 
to the I = m = 2 current multipole: 

r g « -47M^RIaP- 3 s, (9) 
as well as the bulk and shear viscosity damping timescales: 



n 



2.7 x 10 11 Mi.4 J Rr :l P-3T 9 ~ 6 s 



6.7 x 10 7 M~ 



5/4^23/4^2 



10 



(10) 

(ii) 



where we have used the notation Mi. 4 = M /IAMq , Rio = 
R/10 km, P_ 3 = P/l ms and T 9 = T/10 9 K. We have used 
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these timescales because it allows us to make a direct com- 
parison with the earlier work of Owen et al (1998). However, 
it should be noted that this is probably the most optimistic 
scenario with regard to the lifetime of unstable r-modes. We 
do not take into account any damping effects that may arise 
from the interaction of mode and magnetic field (Rezzolla, 
Lamb & Shapiro 2000), or from the presence of exotic parti- 
cles such as hyperons (Jones 2001: Lindblom & Owen 2002). 

2.2 Torques 

The total torque on the star J is taken to be the sum of 
torques due to gravitational radiation, accretion and mag- 
netic dipole radiation. The torque due to gravitational-wave 
emission from the I = m = 2 current multipole is 

j g =3Jtla 2 MR 2 Tg 1 . (12) 

In order to model the accretion torque we need to know the 
quantity and rate of fallback onto the young neutron star. 
Several authors have examined this problem and concluded 
that a significant quantity of matter (up to ~ O.IMq) could 
fall back onto the neutron star during the first few months 
of its life (Colgate 1971; Chevalier 1989; Houck & Chevalier 
1991). Accretion of such a large amount of matter onto a 
neutron star is hypercritical; that is, orders of magnitude 
above the Eddington limit (compare with Bethe, Brown & 
Lee 2000). Such rapid accretion would be advection domi- 
nated, with energy loss occurring via neutrino emission as 
matter is processed at the surface of the neutron star. Ad- 
vection dominated accretion flows have been examined in 
detail by Narayan & Yi (1994,1995). They conclude that 
advection could dominate in cases where the mass accretion 
rate is very high; precisely the conditions likely to prevail 
around a newly-formed neutron star. 

Advection dominated accretion onto a magnetized neu- 
tron star has been examined by Menou et al (1999). Build- 
ing on earlier work (for example Illarionov & Sunyaev 1975; 
Ghosh & Lamb 1978), they model a rotating magnetized 
neutron star surrounded by a magnetically threaded accre- 
tion disk. Accreting matter follows magnetic field lines and 
gives up angular momentum on reaching the surface, exert- 
ing a spin-up torque. The contribution of the magnetically 
threaded disk is however more complex. Interaction between 
the stellar field and the disk results in a positive torque for 
small radii, where the field lines rotate more slowly than the 
local Keplerian speed of the gas. For larger radii, the field 
lines rotate more quickly than the local Keplerian speed, 
resulting in a negative torque. If the spin period becomes 
very short, or the rate of fallback of material onto the mag- 
netosphere drops, the star can enter what is known as the 
propeller regime. Then the rapidly rotating magnetosphere 
exerts a centrifugal barrier that inhibits further accretion. 
Accreting matter is flung away from the star and the star 
experiences a spin-down torque. An accreting magnetized 
star can therefore experience either a spin-up or a spin-down 
torque. 

To proceed we need to determine the likely rate and 
quantity of supernova fallback. Mineshige et al (1997) model 
the supernova fallback problem assuming that the fallback 
material possesses some angular momentum (compare with 
Chevalier (1989) and Brown & Weingartner (1994), who 
model the problem under the assumption of Bondi spherical 



accretion). Mineshige et al find a fallback rate Mf oc t n , 
with 1 < n < 2. We therefore use 



M 



f 



Ct 



-3/2 



(13) 



where C is set by choosing the total mass to be accreted 
M acc between some initial time after the supernova ex- 
plosion and t = oo. The accretion rate is related to the 
fallback rate by 



M 



M f Propeller off, 
Propeller on. 



(14) 



Following Chevalier (1989) and Mineshige et al (1997) we 
choose Macc to be < O.IMq. We choose U = 100 s as this is 
the time after the supernova at which we expect there to be 
a recognisable compact object (Burrows & Lattimer 1986). 
If the propeller effect is not active, we assume that hyper- 
critical accretion continues until the accretion rate falls to 
M < lO~ 4 M0yr -1 . This corresponds to the cut-off point be- 
low which hypercritical accretion is no longer possible. We 
follow Chevalier (1989) and Bethe et al (2000) in assum- 
ing that below this rate photons begin to diffuse out of the 
shocked region and further accretion proceeds at (or below) 
the Eddington rate (~ lO" 8 M yr _1 ). 

If on the other hand the propeller effect becomes ac- 
tive, the situation is less clear. The propeller effect prevents 
accretion and trapped photons may escape as soon as the 
propeller becomes active. But would the pressure of these 
photons reduce the rate of fallback of material on to the mag- 
netosphere? We model two possibilities: i) continued expo- 
nentially decaying fallback, and ii) constant accretion at the 
Eddington rate. For stars with magnetic fields of < 10 13 G 
there is little quantitative difference between the two mod- 
els, however there are noticeable differences for higher fields. 

A typical young neutron star is expected to have a mag- 
netic field B ~ 10 12 G (see Table although a magnetar 
could have a field strength as high as 10 15 G (see for example 
Thompson & Duncan 1996). We assume that the magnetic 
field strength is not affected by accretion or the action of the 
mode, even though it is by no means clear that this would 
be the case in practice. Several authors have suggested that 
the r-mode might lead to differential rotation and hence to 
changes in the magnetic field (see for example Spruit 1999 
and Rezzolla et al 2000). In the earliest stages after the su- 
pernova the accretion rate is likely to be so high that the 
magnetic field of the neutron star is completely confined. 
The pressure of accreting material will be so much greater 
than the magnetic pressure that the magnetospheric radius, 
defined by 



B 2 R b 



2M(2GM) 1 /2 



2/7 



(15) 



with B taken to be the polar magnetic field of the star, is 
smaller than the radius of the neutron star. The magnetic 
field remains confined to the neutron star until the accre- 
tion rate falls sufficiently for the field to begin to exert an 
influence outside the star. 

The accretion torque J a will vary depending on whether 
the magnetic field is confined and whether the propeller ef- 
fect is operational. When the field is completely confined we 
use Narayan and Yi's (1994,1995) results for the accretion 
torque. Their models of advection dominated accretion flows 
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suggest orbital velocities with v 2 



7 u k> 



Vk being the Kepler 



velocity of a particle orbiting the star. Hence 
"2GM1 1/2 



J a = MR 



7R 



(16) 



As the accretion rate falls, r m grows and eventually exceeds 
the neutron star radius 7?. At this stage we switch to the 
equations of Menou et al (1999), as the torque will now 
be affected by the interaction of the magnetic field and the 
fallback material. 



Ja — 2r m i^/ c (^ m ) 



1 - 



0, k {r m ) 



M 
M f 



Propeller off, 
Propeller on, 



(17) 



where fifc(r m ) is the Keplerian angular velocity evaluated 
at the magnetospheric radius. This magnetized accretion 
torque falls as the star spins up, and for f2 > Qk(r m ) the 
torque will in fact be negative. The star is then in the pro- 
peller regime. During this phase matter is prevented from ac- 
creting onto the star and fallback material is instead expelled 
from the system by the rotating magnetosphere. Menou et al 
(1999) review the assumption that accretion ceases and con- 
clude that accretion is indeed minimal when the propeller is 
active. We do not attempt to model the complicated transi- 
tion regime between the two phases in detail. 

Once magnetic field confinement ends we could also in- 
clude a torque due to emission of magnetic dipole radiation, 
J m , where 

JU--^. (18) 

This term was however found to be negligible compared to 
the other torques acting on the young neutron star. Its ef- 
fects will not be discussed further in this paper. 



2.3 Temperature evolution 

Viscosity, the main agent damping the r-mode in our model, 
depends critically on temperature. We include three factors 
in modelling the temperature evolution: modified URCA 
cooling, shear viscosity reheating, and accretion heating. We 
assume the star to be isothermal even though this may not 
be appropriate for a newly born accreting star (Burrows & 
Lattimer 1986). Still, this assumption simplifies the analysis 
considerably and all previous studies of the r-mode instabil- 
ity have assumed a uniform temperature distribution. 

The primary cooling mechanism for a young neutron 
star will be the modified URCA reaction, in which neutrons 
decay via the weak interaction, emitting neutrinos. We do 
not model any rapid cooling effects due to the direct URCA 
reactions in the core of the star (compare with Lattimer et al 
1994 and Page et al 2000). The cooling rate due to the mod- 
ified URCA reaction, i u , is given by Shapiro & Teukolsky 
(1983): 



7.5 x 10 39 Mi/ 4 3 T 9 8 erg s _1 . 



(19) 



The neutron star will be heated by the action of shear viscos- 
ity on the r-mode oscillations. The heating rate due to shear 
viscosity, e s , is given by Andersson & Kokkotas (2001): 



2a 2 Q 2 MR 2 J 



8.3 x 10 3 Vn 2 JM?/ 4 4 fl- 16 / 4 T 9 



erg s 



(20) 



Accretion heating will have two components. The first con- 
tribution arises when accreting matter undergoes nuclear 
burning at the surface of the star. We assume that accret- 
ing matter will release approximately 1.5 MeV of energy per 
nucleon (Brown & Bildsten 1998). If we assume that every 
nucleon reaching the surface is burnt, then energy is liber- 
ated at a rate: 



— x 1.5 MeV 

m B 

4 x 10 51 Mi.4 erg s" 1 



(21) 



where tub is the mass of a baryon. 

A second contribution arises because the flow is as- 
sumed to be advection dominated. Matter falling in towards 
a star will liberate ~ GM/R of potential energy per unit 
mass. In a non-advection dominated flow most of this en- 
ergy would be dissipated as heat in the accretion disk. In an 
advection dominated flow, by contrast, this energy is carried 
in with the flow of matter. We assume that the advected po- 
tential energy is used to generate neutrinos near the surface 
of the star, and that these neutrinos are radiated isotropi- 
cally from the point of generation. We therefore assume that 
half of the neutrinos escape without interacting with the star 
and that the remaining half are radiated into the star, where 
they scatter and interact with the stellar material. The mean 
free path A of inelastic scattering of neutrinos with electrons 
in the neutron star matter is 



A~2x 10 



where p„ 



pnuc 


7/6 


rO.l MeV] 


9 




L E„ \ 



-,/2 



(22) 



2.8 x 10 14 g cm 3 and p is the average density 



of the star. For an n — 1 polytrope we have 



P 

pnuc 



0.7 



M 

Mo 



(23) 



We assume that the neutrino energy is E v — 1 MeV (Shapiro 
& Teukolsky 1983). Assuming fully efficient scattering, and 
noting that A > R, the heating rate due to advected poten- 
tial energy may be estimated as 



R GMM 
1 R 

= 8 x 10 51 M^ 6 Mia erg s" 



(24) 



We estimate the heat capacity C„ by assuming the thermal 
energy of the neutron star to reside almost exclusively in 
degenerate neutrons (Shapiro & Teukolsky 1983). Neglecting 
interactions, the heat capacity of such a system is given by 



C v = 1.6 x 10 39 M^ 4 3 T 9 erg K" 



(25) 



Equation (|25|) shows that C„ oc T, hence the equation of 
thermal balance of the star is 



-C V T 



S u + in + Sh 



(26) 



Manipulation of equation (|26|) yields an expression for the 
temperature evolution of the star. We initiate our evolutions 
at T — 10 11 K and terminate them at T = 10 9 K, the 
temperature at which we expect a crust to form (Haensel 
1997) and suppress the r-mode instability in an accreting 
star. 
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3 OBSERVATIONAL CONSTRAINTS 

Before we proceed to the discussion of the results obtained 
from our simple model it is useful to consider the constraints 
provided by observations of young neutron stars. It is stan- 
dard practice to assume that a pulsar spins down purely as 
a result of magnetic dipole radiation. Given the observed 
period and the spin-down rate we can then estimate the 
magnetic field strength from 

B = 3.2 x 10 19 (PP) 1/2 G. (27) 

Similarly, one can infer the characteristic age t c of the pulsar 
from the observed data: 
P_ 
2P~ 

Here it is assumed that the star was born spinning much 
faster that its current rate. The use of this estimate is sup- 
ported by the fact that, in the case of the Crab pulsar, the 
age calculated in this way accords well with the known age. 

If observations also provide the braking index, n, and 
an estimate of pulsar age, tsNR, the initial spin period of 
the pulsar Po is inferred from 

l/(n-l) 



tc = 



(28) 



Po = P 



(n- 



P 



where n is given by 



PP 



2. 



(29) 



(30) 



In most cases the largest uncertainty in estimating the 
initial spin period is associated with the age of the pulsar, 
tsNR- Ideally, one would like to make an association with a 
known historical supernova. This would then, as in the case 
of the Crab, provide a precise age estimate and a reason- 
ably reliable initial spin period. Unfortunately, this kind of 
data is the exception rather than the norm. Typically, the 
independent age is estimated by measuring the size of the 
supernova remnant and comparing it with theoretical mod- 
els for the expansion rate. This procedure is complicated by 
many factors and the results are associated with large error 
bars. 

A sample of data for young pulsars that are claimed 
to be associated with known supernova remnants is given in 
Table [l]. The reliability of the inferred initial spin periods can 
be assessed by comparing the characteristic age of the pulsar 
to the estimated age of the supernova remnant. For most of 
the cases listed in Table [j], these estimates are consistent. In 
these cases the association between supernova remnant and 
pulsar seems correct and the spin of these pulsars appears to 
have evolved primarily due to the action of magnetic dipole 
radiation. 

The data for young pulsars clearly support the notion 
that neutron stars are not formed spinning near the Ke- 
pler limit. This would be consistent with both the r-mode 
instability scenario (where a rapidly spinning, newly born 
neutron star is spun down to a period of say 10 ms in a 
few months) and the possibility that magnetic core-envelope 
coupling in the supernova progenitor leads to most neutron 
stars being born rotating slowly (Spruit & Phinney 1998). 

There is ongoing discussion in the literature as to 
whether the (obviously simplistic) magnetic dipole model 
adequately describes the spin evolution of young pulsars. 



Concerns stem from discrepancies between various age es- 
timates for some young pulsars. Calculations would also be 
complicated by magnetic field decay (Colpi et al 2000) or 
growth (Gaensler & Frail 2000). Another problem is that all 
measurements of the braking index n are lower than the 
value of 3 predicted by the magnetic dipole model. Sev- 
eral authors have suggested that the propeller effect, driven 
by ongoing low rate accretion of supernova remnant ma- 
terial, could explain the discrepancies (Marsden, Lingenfel- 
ter & Rothschild 2001; Menou, Perna & Hernquist 2001). 
The same mechanism might also explain the properties of 
the Anomalous X-Ray Pulsars (AXPs) (Perna, Hernquist & 
Narayan 2000; Chatterjee, Hernquist & Narayan 2000; Al- 
par 2001; Menou, Perna & Hernquist 2001), without having 
to invoke large magnetic field strengths (Thompson & Dun- 
can 1996). See Rothschild, Lingenfelter & Marsden (2001) 
for a review of AXP properties and the competing accretion- 
driven and magnetar models. 

The relation of these propeller /fallback models to the 
one used in this paper is simple. Because we are interested 
in gravitational wave emission, we model only the very ear- 
liest stages of fallback and terminate our evolutions at crust 
formation. At crust formation (t ~ 1 year) there is still ma- 
terial in the fallback disk, although the fallback rate has 
dropped significantly. The papers listed above focus more on 
the subsequent long period of sub-Eddington accretion. Al- 
though they do discuss the early super-Eddington accretion 
phase (see for example Chatterjee, Hernquist & Narayan 
2000), they do not model this phase in detail and neglect 
the effects of gravitational wave emission completely. 



4 RESULTS FROM THE FALLBACK MODEL 

We are now prepared to discuss the results obtained from 
our phenomenological r-mode instability/supernova fallback 
mode, and compare the predicted spin rates to the data in 
Table [lj. 

If we consider the data in Table [j], and the fact that the 
r-modes are certainly stable for periods longer than (say) 
20 ms (see Andersson & Kokkotas 2001), it would seem as 
if the instability may only be relevant for a small sample 
of nascent neutron stars. This would be in agreement with 
Spruit & Phinney 's (1998) model, which suggests that most 
neutron stars are born rotating slowly. Our study will show, 
however, that supernova fallback may change this picture 
considerably, with accretion induced spin-up driving the star 
into the window where the r-mode is likely to be unsta- 
ble. In other words, significant supernova fallback enhances 
the probability that a young neutron star evolves through a 
phase where the r-mode instability is active. 

Let us consider the effect of supernova fallback on a 
neutron star born spinning at a reasonably small fraction of 
the Kepler rate. Taking the birth spin period as 10 ms, we 
evolve the star according to our phenomenological model. 
Typical results for the resultant spin evolution are shown in 
Figure jjj. The data shown in the figure suggest the following 
picture. After an initial spin-up phase, the r-mode becomes 
unstable and begins to grow. For a typical young neutron 
star with magnetic field ~ 10 12 G the r-mode instability pro- 
vides the main spin-down mechanism. For larger magnetic 
fields the propeller effect becomes active at an earlier stage 
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Figure 1. The solid lines show the evolution of the spin period for 
newly-born accreting neutron stars with an unstable r-mode and 
different magnetic field strengths [Key: Hollow circles - propeller 
switches on; Upwards triangles - r-mode saturates; Downwards 
triangles - r-mode becomes stable]. Note that the r-mode never 
reaches the saturation amplitude for the star with the highest 
magnetic field. The dashed lines, given for comparison for the 
two extreme magnetic field strengths, show the spin evolution if 
an unstable r-mode is not present. 
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Figure 2. Spin rates predicted by the fallback model at crust 
formation (within 1-10 years of the supernova) compared to val- 
ues inferred for young pulsars using the standard magnetic dipole 
spin-down model (black squares). The solid line is for a model 
where fallback rate continues to decay exponentially when the 
propeller switches on; the dashed line for one where the fallback 
rate then drops to the Eddington limit. [Note that the obser- 
vational data is associated with large error bars which are not 
included in the figure.] 



of the evolution and is the main spin-down mechanism for 
highly magnetized stars. Nevertheless, the r-mode is impor- 
tant in even the most highly magnetized stars. Switching 
off the r-mode leads to a marked reduction in spin-down, 
despite the fact that the mode never actually saturates in 
the star with the largest magnetic field. Note also that our 
model predicts that supernova fallback will delay the ap- 
pearance of the r-mode signal compared to the case of an 
isolated star, for which the r-mode is predicted to saturate 
after ~ 10 minutes (Owen et al 1998). This can be seen from 
equation (^): the high initial accretion rate overwhelms all 
of the other terms in the expression and a will not grow 
until M has fallen. 

In our model we have assumed that the magnetic field 
does not influence Mj (although it does affect M via the 
propeller effect). Recent work by Igumenshchev & Narayan 
(2002) on spherical fallback, however, suggests that the pres- 
ence of a magnetic field may reduce Mf by perhaps a factor 
of 10. If we assume such a reduced Mf from the time at 
which magnetic field suppression ends, we find little impact 
on spin rate or gravitational wave emission for a 10 12 G star. 
The effect on a 10 14 G star is however more pronounced. Ini- 
tial spin-up and subsequent propeller-driven spin-down are 
much reduced, as is the emitted gravitational wave ampli- 
tude. Note however that we only examine the first ~ 1 year 
post-supernova and that the effects of slower fallback may 
be more important at later times. 

We have terminated the evolutions once the star reaches 
a temperature of 10 9 K. At this point one would expect 
the crust to form, which is predicted to increase dramati- 
cally viscous damping of the r-mode (Bildsten & Ushomirsky 
2000; Lindblom et al 2000). In Figure | we compare the spin 
period at the end of the our evolutions to the values of Pn in- 
ferred for young pulsars using the magnetic dipole model (as 
listed in Table hi). This figure should, of course, be considered 
with caution since the inferred spin periods are associated 
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Figure 3. The combined action of the r-mode instability and 
fallback accretion causes the spin rates of stars with birth spin 
rates in the range 1 ms to 1 s to converge to ~ 10 ms within ~ 1 
year. The data shown are for a star with B = 10 12 G. 

with significant error bars. Still, our model clearly captures 
the main trend of the data, namely that the initial spin pe- 
riod inferred from the magnetic dipole model increases with 
magnetic field. Given the many uncertainties in both the 
theoretical model and the observational data the agreement 
between the results is rather good. 

The fallback model suggests that the outcome is insensi- 
tive to the spin period at the beginning of the evolution (the 
true "birth" spin rate). The combined action of supernova 
fallback and the r-mode instability causes the spin rates of 
stars with very different birth spins to converge within ~ 1 
year. This is shown in Figure [| 

This result is notable given the uncertainty over the 
spin rates at which neutron stars are formed. It may, in fact, 



The spin evolution of nascent neutron stars 7 



20 -, 




-) 

100 1000 10000 100000 1000000 1E7 1E8 1E9 

Time (s) 



Figure 4. The effect on spin evolution of changing a B for a star 
with B = 10 12 G. 

suggest that the true "birth" spin is not the key factor that 
determines the spin period at (say) one year post-supernova. 
This observation has implications for determining whether 
or not a young neutron star undergoes a phase where the 
r-modes are unstable. This will in turn affect the relevance 
of the r-modes as a gravitational wave source. Of particular 
interest is the fact that our results suggest that fallback 
accretion may lead to even a very slowly spinning neutron 
star being spun up sufficiently for the r-modes to become 
unstable after ~ 10 min. However, we conclude that the 
magnetic field is far more important in determining the spin 
rate at one year post-supernova than both the birth spin 
rate and the possible r-mode spin-down phase. 

The key parameter in determining the role of the r- 
mode instability is the saturation amplitude a s . As has 
been demonstrated by recent hydrodynamical simulations 
(Stergioulas & Font 2001, Lindblom, Tohline & Vallisneri 
2001), the r-modes may not saturate until the amplitude 
has reached values of order unity. The effect of varying a a is 
shown in Figure ^[ These results are exactly what one would 
expect: the r-mode needs to be able to grow to a large ampli- 
tude if it is to counteract the strong accretion torque during 
the fallback phase. 

One might expect fallback accretion to lead to strong 
heating, and that this would have a significant effect on the 
duration of the r-mode gravitational-wave signal. We anal- 
yse the role of accretion heating in Figure ]A We find that 
even though rapid accretion increases the temperature by 
roughly 10 percent throughout the initial evolution, it does 
not lengthen the r-mode instability phase. Far more impor- 
tant in this respect is reheating due to the action of shear 
viscosity on the r-mode, which prolongs the r-mode insta- 
bility phase for both an accreting and a non-accreting star. 
The amount of viscous heating is essentially the same in our 
current model as it would be in the case of an isolated young 
neutron star. 

In order to assess the detectability of the gravitational 
waves due to the r-modes, we use the prescription of Owen 
et al (1998). The true strain of the gravitational waves, h, 
is given by 

h(t) = 7.54 x W' 23 aJfl 3 M 1A Rl 15 ^ pC . (31) 
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Figure 5. The effect of shear viscosity reheating and accretion 
heating on the temperature evolution of a young neutron star. 



We consider sources in the Virgo Cluster, for which the dis- 
tance to the source, D ~ 15 Mpc. Results from our model 
are similar to results from previous studies (such as Owen 
et al 1998). They show that the unstable r-modes may lead 
to a gravitational wave signal that could be detectable by 
a second generation detector like LIGO II. The main differ- 
ence between our results and those of Owen et al (1998) is 
the duration of the signal. We predict that, if continued ac- 
cretion were to take place and prevent crust formation until 
the core temperature has fallen to 10 9 K, the signal may 
last for several years. We find that the peak gravitational- 
wave signal amplitude remains almost unchanged as we vary 
the saturation amplitude by two orders of magnitude, even 
though this has a noticeable effect on final spin period. Note 
however that these conclusions are only relevant for neutron 
stars with fields < 10 13 G. It is unlikely that we would ob- 
serve a signal from the r-modes of neutron stars with larger 
magnetic fields. This conclusion agrees well with studies of 
Rezzolla et al (2000) and Ho & Lai (2000). In our case the re- 
duction in signal with increasing magnetic field is due mainly 
to increased propeller effect spin-down. This is because the 
gravitational wave amplitude depends so strongly on spin 
rate (see equation (|3l|)) . It should of course be pointed out 
that we have not taken account of any possible interaction 
between the r-mode and the magnetic field, which may fur- 
ther decrease the relevance of the r-modes in magnetized 
stars. 

Finally, we have considered the possibility that the fall- 
back accretion torque might be able to spin the star up 
to (and beyond) the mass-shedding limit (for the canonical 
n — 1 polytropic neutron star used in our evolutions this cor- 
responds to P « 0.8 ms). Depending on the overall stiffness 
of the supranuclear equation of state, neutron stars spinning 
near the mass-shedding limit may be subject to the dynam- 
ical bar-mode instability (Shapiro & Teukolsky 1983). This 
would lead to the formation of a highly non-axisymmetric 
bar-like configuration which is expected to be a strong source 
of gravitational waves with frequencies of order kHz (Houser, 
Centrella & Smith 1994). Whether a young neutron star be- 
comes bar-mode unstable is therefore of great interest to 
gravitational wave theorists. Interestingly, our evolutions in- 
dicate that the r-mode instability is remarkably efficient in 
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preventing the star from reaching the mass-shedding limit 
during the fallback phase. Even stars born with periods as 
short as 0.85 ms were unable to reach the Kepler thresh- 
old. The r-mode grows rapidly and is able to spin down the 
star, even in the initial stages where the accretion torque is 
the highest. This suggests that r-mode induced spin-down 
might counteract accretion spin-up to prevent all but the 
most rapidly rotating young neutron stars from becoming 
bar-mode unstable. For comparison we also investigated the 
case where there was no unstable r-mode. In this case the 
propeller effect switches on far earlier and can also act to 
prevent the star from going bar-mode unstable. Consider 
the B = 10 12 G star shown in Figure [l], which is born spin- 
ning at P = 10 ms. If the unstable r-mode is present, the 
propeller effect will not switch on until ~ 200 days. If the 
r-mode is absent, however, the propeller will come on after 
only ~ 20 days. For a B = 10 12 G star born spinning at 
P — 1 ms, the difference is even more pronounced. In this 
case the propeller will switch on after ~ 3 days when there is 
no r-mode, as compared with ~ 200 days when the unstable 
r-mode is present. 



5 CONCLUDING REMARKS 

We commenced this work with two prime objectives: to un- 
derstand the impact of supernova fallback accretion on the 
unstable r-modes of neutron stars and any associated gravi- 
tational wave signal, and to investigate the combined influ- 
ence of accretion and r-mode instability on the spin rates of 
young neutron stars. Our model, simple as it is, has offered 
insight into both questions. Our results suggest that newly 
formed neutron stars accreting fallback material will indeed 
experience a phase of perhaps several years during which 
the r-mode is unstable. This would have positive ramifica- 
tions for detection of gravitational waves from the unstable 
r-modes of neutron stars. Moreover, the model suggests that 
this is true even for stars that are born with low spin rates. 
The effect of magnetic field is however crucial. While stars 
with canonical (< 10 13 G) field strengths are predicted to go 
r-mode unstable, the strong propeller effect associated with 
magnetars appears to preclude the r-mode from growing to 
the saturation level. The external torque simply spins the 
star down fast enough that the r-mode is not given time to 
grow. It is thus unlikely that we will see a strong r-mode 
gravitational wave signal from a highly magnetized star. 
The likelihood of obtaining a gravitational wave signal from 
an unstable bar-mode is also low in our model. Spin-down 
caused by emission of radiation from the r-modes is able to 
counteract even the strongest accretion torque and prevent 
all but the stars with the most rapid initial spins (and obvi- 
ously those born spinning faster than the Kepler limit) from 
going dynamically unstable to the bar-mode. The propeller 
effect has a similar preventive effect for the stars with the 
highest magnetic fields. 

With regard to our second question, the model's pre- 
dictions are consistent with spin rates inferred from obser- 
vations of young pulsars. A canonical neutron star is spun 
down within a few years by the combined action of fall- 
back and r-mode to a period similar to that inferred for the 
Crab pulsar at its birth. Neutron stars with stronger mag- 
netic fields, for which the propeller effect is the dominant 



spin-down mechanism, achieve far longer periods in the first 
year post-supernova. Even for such highly magnetized stars, 
however, the effect of an unstable r-mode on spin rate is 
significant. 
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Table 1. Young pulsars (PSR) that have been identified with known supernova remnants (SNR). Where the braking index is not provided 
by observation, it is assumed to be similar to that of the Crab (n = 2.5). The lower part of the table illustrates the difficulty of estimating 
tSNR • 

Notes: (a) Measurement technique for n is somewhat unorthodox, (b) Age estimate made using proper motion of pulsar, (c) Using 
tsNR = 900 years - for larger isjvij it is not possible to estimate Po using equation (29). (d) As tg^R = t c it is not possible to estimate 
a value of Po using equation (29) - we can only say that Po is small, (e) Age of supernova remnant too uncertain to make a reliable 
estimate of Po. (f) Estimated age of remnant much larger than t c making estimate of Po unfeasible. 

Most references from Kaspi (2000), additional data from: B0833-45 (Lyne et al 1996), B0540-69 (Zhang et al 2001, Reynolds 1985), PSR 
1951+32 (Foster et al 1994, Migliazzo et al 2002), J1811-1926 (Torii et al 1999), J1846-0258 (Gotthclf et al 2000), J205+6449 (Murray 
et al 2001), J1124-5916 (Camilo et al 2002), B1853+01 (Cox et al 1999), B1509-58 (Kaspi ct al 1994), J1119-6127 (Crawford et al 2001), 
J1016-5857 (Camilo ct al 2001), B1757-24 (Gaenslcr & Frail 2000). 



PSR 


P (ms) 


n 


SNR 


SNR age (kyr) 


tc (kyr) 


Po (ms) 


B (10 12 G) 


B0531+21 


33 


2.52 


Crab 


0.95 (SN1054) 


1.2 


18 


3.8 


B0833-45 


89 


1.4 (a) 


Vela 


9-30 


11 


13-57 


3.4 


B0540-69 


50 


1.81 


N158A 


0.8-1.1 


1.7 


34-39 


5.0 


PSR 1951+32 


39.5 




CTB 80 


38-74 (b) 


107 


21-33 


0.49 


J1811-1926 


65 




Gil. 2-0.3 


1.7 (SN386) 


24 


63 


1.7 


J1846-0258 


324 




Kcs 75 


0.9-4.3 


0.7 


315 (c) 


49 


J205+6449 


66 




3C58 


0.82 (SN1181) 


5.4 


60 


3.6 


J1124-5916 


135 




G292.0+1.8 


1.7 


2.9 


90 


10 


B1853+01 


267 




W44 


20 (d) 


20 


« P 


7.6 


B1509-58 


151 


2.83 


MSH15-52 


(e) 


1.6 


(e) 


15 


J1119-6127 


408 


2.91 


G292-05 


(e) 


1.6 


(e) 


41 


J1016-5857 


107 




G284.3-1.8 


~ 10 (c) 


21 


(e) 


3 


B1757-24 


125 




G5.4-1.2 


39-170 (f) 


15.5 


(f) 


4 



